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A method for the estimating of the residual service life of an apparatus 



10 The invention relates to a method for the estimating of the residual service 
life of an apparatus which is subject to wear during operation. The inven- 
tion further relates to the use of such a method for service planning. 

For many apparatuses, such as turbines or jet engines of aeroplanes, the 
1 5 wear of or at important components is the main reason for the limited 
service life, with "service life" usually meaning the period of time or the 
operating period between two standard audits. Wear can manifest itself by 
changes in the mechanical properties of components, which is due, for 
example, to friction, heat or fatigue. These changes at the components 
) 20 also bring about a changed behaviour of the whole apparatus with un- 
changed operating conditions. It is therefore necessary to carry out audits 
or service work on the apparatus at regular intervals. 

Methods known today for the planning of such service work are usually 
25 based on purely time-defined service intervals, that is a fixed time interval 
is set for the service life between two audits. For safety reasons, very 
conservative approximations are assumed for the service life already used 
up. It is clear that such purely time-supported methods do not result in 
an optimum utilisation of the possibilities of the apparatus, because the 
30 wear which actually occurred, which is also influenced by the specific 
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operating conditions and environmental conditions, is not taken into 
account. To increase the efficiency of utilisation of such apparatuses, for 
example aeroplane engines, it is therefore desirable to have a more realis- 
tic estimate available for the residual service life, in which the actual 
5 operating conditions and environmental conditions are also taken into 
account in order thus to maike possible a more efficient and more eco- 
nomic service planning. 

J The present invention is directed to this object. A method should therefore 

10 be provided for the estimating of the residual service life of an apparatus 
in which the specific operating conditions, and thus the service life of the 
apparatus effectively, i.e. really, used up, are taken into account. 

The method satisfying this object is characterised by the features of the 
15 independent method claim. 

In accordance with the invention, a method is therefore provided for the 
estimating of the residual service life of an apparatus which is subject to 
wear during operation, having the following steps: 
) 20 a) for at least one characteristic parameter which is sensitive to the wear, 
a relationship is determined to a time parameter which is representative 
for the operating period; 

b) a limit value is fixed for the characteristic parameters which gives the 
maximum permitted wear; 
25 c) a code field is established which gives a relationship between the char- 
acteristic parameter, the time parameter and the wear; 
d) actual values are determined for the characteristic parameter in de- 
pendence on the time parameter with the aid of data obtained by a techni- 
cal measurement; 
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e) the instantaneously present wear is determined from the actual values 
with reference in each case to the code field; 

f) starting from the instantaneous actual value of the characteristic pa- 
rameter, a determination is made by means of extrapolation to the limit 

5 value of the end value of the time parameter for which the maximum 
permitted wear is reached; 

g) the residual service life is estimated by a comparison of this end value 
with the value for the time parameter which belongs to the instantane- 

^ ously present wear. 

10 

In the method in accordance with the invention, the wear is quantified 
with reference to at least one characteristic parameter which is sensitive 
with respect to the wear. A code field is established for the apparatus 
which describes the relationship between the characteristic parameter, the 
15 time parameter and the wear. In the case of only one characteristic pa- 
rameter, this code field can be represented as an area in a three- 
dimensional space which space is set up by the characteristic parameter, 
the time parameter and the wear. 

) 20 Data are detected at the apparatus in a technical measurement from 

which an actual value can be determined for the parameter in dependence 
on time. A determination can then be made with reference to the code field 
as to how far the wear has progressed in a quantitative manner, for exam- 
ple in per cent. The residual service life can then be estimated by extrapo- 
25 lation to the maximum permitted limit value for the characteristic parame- 
ter. 

The method in accordance with the invention thus takes the service life 
into account which has been actually and effectively used up in order to 
30 estimate the residual service life starting from this. Not only the time 
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parameter is thus taken into account, but in addition, the operation con- 
ditions, and optionally the environmental conditions are taken into ac- 
count under which the apparatus was operated up to the present point in 
time. This operation dependent estimate of the residual service life, which 
5 takes the history of the apparatus into account, makes a substantially 
more efficient utilisation of the apparatus possible, because services are 
only carried out when they are actually necessary. The reliable prediction 
of the wear development thus makes a condition-based service planning 
possible. 

10 

The residual service life is preferably determined for different stages in the 
life cycle, that is, as the life is progressively used up - measured by the 
time parameter - the residual service life is repeatedly re-estimated. 

15 In a particularly preferred embodiment, the code field is established with 
the aid of a-priori knowledge of the wear behaviour. Since apparatuses 
such as aeroplane engines (jet engines) represent very complex systems, it 
is as a rule only possible - if at all - with a fairly great effort to carry out a 
sufficiently precise physical or deterministic modelling of the apparatus. It 

20 is therefore preferred to use a-priori know-how for the establishing of the 
code field. The experience, observations or also measurements which have 
been collected at the same or similar apparatuses are used to describe the 
qualitative and/ or quantitative behaviour of the characteristic parameter 
as the wear progresses. The establishing of the code field is made substan- 

25 tially simpler by the use of such a-priori knowledge. Furthermore, the use 
of the a-priori knowledge normally results in the code field describing the 
apparatus better or more exactly. 
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The a-priori knowledge preferably includes the qualitative and/ or quanti- 
tative course of wear curves which give the relationship between the char- 
acteristic parameter and the time parameter. 



5 Since the a-priori know-how usually includes specifications available in 
verbal form, the code field is particularly preferably established by means 
of a linguistic fuzzy model. A purely qualitative model, for example, can 
thereby be generated to determine the wear, said model then being ade- 
quately optimised with respect to its quantitative properties on the basis 
10 of measured life cycles. 

It has also proved advantageous in practice for the code field to be modi- 
fied with reference to measurement data or on the basis of plausibility 
observations. Such plausibility observations, for example, have proven to 
15 be very useful to model certain marginal regions of the code field which 
correspond to states which the real apparatus seldom reaches. 



It is furthermore preferred for the data obtained by a technical measure- 
ment each to be subjected to a filtering or to an averaging for the determi- 
) 20 nation of the actual values for the characteristic parameter. The data 

obtained by a technical measurement are frequently overlaid by a noise or 
another interfering value such that their direct use, in particular with 
fuzzy models, does not allow any sound statement with respect to the 
degree of wear. 

25 

For the determination of the actual values, a model is established, pref- 
erably with the aid of a plurality of sets of data obtained by a technical 
measurement, with which an actual value is determined for the character- 
istic parameter. 
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The metJiod in accordance with the invention is in particular suitable for 
the estimating of the residual service life of an engine, in particular of an 
aeroplane engine. 

5 The method in accordance with the invention is particularly well suited for 
the service planning, in particular of an aeroplane or of a plurality of 
aeroplanes (fleet management) . 

Further advantageous measures and preferred aspects of the invention 
10 result from the dependent claims. 

The invention will be described in more detail in the following with refer- 
ence to an embodiment and to the drawing. There are shown in the sche- 
matic drawing: 
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Fig. 1: 



typical wear curves for an aeroplane engine; 



Fig. 2: 



membership function of a fuz2y set; 



20 Fig. 3: 



a code field in a three-dimensional space which is set up by 
the A axis, the T axis and the V axis; 



Fig. 4: 



a projection of historic data into the plane which is set up by 
the A axis and the V axis; 
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Fig. 5: 



a presentation of the A-T plane for the illustration of an ex- 
trapolation in an embodiment of a method in accordance with 
the invention; 
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Fig. 6: 



as Fig. 4, but additionally with the projection of an extrapola- 
tion; and 



Fig. 7: 



a plurality of representations, in each case as in Fig. 6; how- 
ever, at different times in the life cycle, for the illustration of 
the respective estimate of the residual service life in accor- 
dance with the embodiment of the method in accordance with 
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the invention. 



) 

10 The method in accordance with the invention is explained in the following 
with an exemplary character with reference to a jet engine of an aeroplane. 
The engine stands as a representative example for an apparatus which is 
subject to wear during operation and for which an estimate of the residual 
service life should be made. The invention is naturally not restricted to 

15 this application, but is also suitable in basically the same manner for 



) 20 The term "service life" means the interval between two regular audits or 
services, that is the degree of wear is evaluated as "new" directly after the 
service. 

The change in the gas discharge temperature of the engine is selected as 
25 the chsiracteristic parameter which is sensitive with reference to the wear 
of the aeroplane engine. The temperature change is designated by T in the 
following. 



other apparatuses such as land based turbines, flow machines or other 
mechanical systems which are subject to wear during operation and 



therefore have to be serviced. 
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The number of flights is selected as the time parameter which is represen- 
tative for the operating time or the service life and is designated by A in 
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the following. Other time parameters would naturally also be suitable, for 
example the number of operating hours. 

The wear V is described by a number from the interval from zero to one, 
5 where V = 0 means "new", no wear has yet occurred, and V = 1 means 
"used", that is the permitted wear limit has been reached. A wear of V = 
0.6 therefore means, for example, that the degree of wecir amounts to 60% 
of the maximum permitted wear. 

10 A temperature increase of 40 K is chosen as the limit value G for the 

temperature change T which gives the reaching of the maximum permitted 
wear. This value is based on experience which is brought in as a-priori 
know-how. 

15 Fig, 1 shows a plurality of typical wear curves K1-K4 which each give a 

possible relationship between the temperature change T as the character- 
istic parameter and the number A of flights as the time parameter. Such 
wear curves represent a-priori knowledge which is used in the present 
embodiment for the establishing of the code field. The wear curves are 

20 based on experience or data which are obtained by a technical measure- 
ment and have been made on the same or similar engines. 

Furthermore, the limit value G = 40 K is drawn in for the temperature 
increase as well as an interval L which gives a typical life cycle for such an 
25 engine. 

It can be recognised that the curves Kl - K4 admittedly all look different, 
but show the same behaviour qualitatively. At the start of the life cycle, 
which begins at A = 0 and T = 0, a steep climb of the wear curves can first 
30 be recognised; a plateau follows at which the temperature increase hardly 
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changes; and towards the end of the life cycle a steeper climb in the tem- 
perature increase can again be observed. All wear curves end at T = 40K, 
that is on reaching the limit value G. 

5 A possibility of quantifying the wear consists of the fact of norming the 
path integral over the whole wear curve to one. Then the path length for 
each point on the wear curve defines exactly the corresponding wear. If, 
for example, Al flights have been covered, the quantified wear results from 
the path integral over the wear curve from 0 to Al. 

10 

The problem is, however, that for a given engine the wear curve obtained 
by a technical measurement is only known at the end of its lifetime. Ex- 
pressed differently: it is not known for a new, i.e. freshly serviced engine, 
along which of the many different wear curves it will move during its life 
15 cycle. 

Expressed in somewhat simplified terms, in the embodiment of the 
method in accordance with the invention described here, the estimate of 
the residual service life is improved in the course of the progressive using 
20 up of the service life in that the prediction is constantly (that is, for exam- 
ple, after every flight) matched to the correct wear curve. 

From the viewpoint of an input/ output model, the characteristic parame- 
ter and the time parameter which measures the service life used up to 
25 date are to be used as inputs and the wear as the output as the required 
and adequate information for the estimate of the residual service life. 

In the method in accordance with the invention, a code field KF (see Fig. 3 
is first established which gives a relationship between the characteristic 
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parameter (here tiie temperature increase T), the time parameter (here the 
number of flights A) and the wear V. 

This code field is preferably established with the aid of a-priori know-how, 
5 that is for example using wear curves Kl - K4, as are shown in Fig. 1. 

A linguistic fuzzy model (MAMDANI model) is particularly preferably used 
for the conversion of this a-priori knowledge. Since such fuzzy models or 
^ the fuzzy logic per se are sufficiently known to the person skilled in the 

10 art, it will only be explained briefly here how the code field KF is estab- 
lished with the aid of a linguistic fuzzy model in the embodiment of the 
method in accordance with the invention. 

The following information is available as a-priori know-how for the linguis- 
15 tic fuzzy model in the embodiment described here: knowledge of the quali- 
tative and quantitative course of wear curves (see Fig. 1); additionsd de- 
mands on the prediction in marginal areas. 

The temperature increase T and the service life used up to date, measured 
) 20 by the number A of flights, serve as input vedues, the wear V as the output 
value. Each of the input and output values is characterised by a fuzzy set. 
As an example, Fig. 2 shows the membership function of the fuzzy set for 
the input value temperature increase T. The linguistic values "small", 
"medium" and "large" are provided for the linguistic variable T. 

25 

The linguistic values "few", many", "limit" are provided for the input value 
A "number of flights" and the linguistic values "new", "as new", "used", 
"used up" are provided for the output value "wear" V. 
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Subsequently, the rules are defined using the fuzzy set; in the specific 
case, these includes the following four rules 

IF (T is small) AND (A is few) then (V is new) 
5 IF (T is small) AND (A is many) then (V is as new) 
IF (T is medium) AND (A is many) then (V is used) 
IF (T is large) AND (A is limit) then (V is used up). 

A code field KF can then be established from this information as is shown 
10 as a mesh in Fig. 3. It is naturally possible, and normally also usual, for 
initially a first "draft" to be established for the code field KF and for this 
then to be fine tuned. Such fine tuning operations can take place, for 
example, by trial and error, by taking measurement data into account, by 
manual remodelling on the basis of plausibility observations (e.g. that the 
15 margin of the code field lies on the line with T = 40 K and V = 1) or by 
other optimisation or calibration methods. The code field KF resulting 
from this for the embodiment described here is shown in Fig. 3. 

With the aid of this code field KF, which describes the relationship be- 
) 20 tween the temperature increase T, the number of flights A and the wear V, 
the estimate of the residual service life is now carried out during the life 
cycle of the engine, which will be described in the following. 

During the life cycle of the engine, technical measurement data of the 
25 engine are detected - for example on or after each flight - from which 

actual values are then determined for the characteristic parameter (here 
the temperature increase T) in dependence on the time parameter (here 
the number A of flights). It is naturally also possible, and on occasion also 
advantageous, to use the data detected directly in a technical measure- 
30 ment as the actual values depending on the specific application; in the 
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present embodiment, however, it has proven to be advantageous to sub- 
ject the data obtained by a technical measurement to a filtering or an 
averaging. The reason for this is that the change in the gas discharge 
temperature T is overlaid by a strong noise due to the high absolute value 
5 of the gas discharge temperature. A typical noise amplitude can amount to 
half or even more of the data signals obtained by a technical measure- 
ment. 



Basically, a very large number of processes and methods known per se are 
10 naturally suitable to subject the data obtained by a technical measure- 
ment to an averaging or to a filtering. A possibility is explained with exem- 
plary character which uses a method which is described in EP-A- 895 197 
(P.6822). 

15 The change in the gas discharge temperature T is obtained by a technical 
measurement at regular intervals, that is on every flight, for example. This 
results in a data set of the form [Ti ; Ai], with i = 1, n, where ti is the 
temperature change for the flight with the number Ai. Part of this data 
set, for example the data with i = 1, 2, 20; the data with i = 21 22, 

20 40; etc., is now in each case used to establish a model. This results in a 
plurality of models. Each of these models is then evaluated for the same 
state, the so-called nominal state in order to determine in this manner the 
actual values for the characteristic parameter. Reference is made to EP-A- 
0 895 197 with respect to further details. 
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An averaging of the data obtained by a technical measurement takes place 
by this measure which substantially reduces the noise amplitude. In this 
manner, actual values can be determined for the temperature increase T 
in dependence on the number A of flights. 
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The historic measurement data or the actual values for the temperature 
increase T up to the present are available for the engine at the time at 
which the estimate of the residual service life should take place. A point in 
the plane set up by the A axis and the T axis (see Fig. 3) corresponds to 
5 each value padr (Tk, Ak) from an actual value for the temperature increase 
and the associated service life, measured by the number of flights Ak. For 
reasons of better clarity, only one such point is drawn as a cross in Fig. 3 
with the coordinates (Tk, Ak) in the T-A plane. This point is now projected 
onto the code field KF, that is - in accordance with the representation - 
10 upwardly onto the code field represented by the mesh KF. Then the wear 
Vk associated with (Tk, Ak) can be read off directly. 

As already mentioned, all value pairs (Tk, Ak) up to the present are known. 
If each value pair is imaged onto the code field KF and if the result is 

15 projected onto that plane which is set up by the A axis (number of flights) 
and the V axis (wear), the curve shown in Fig. 4 results, for example. It 
describes the quantified wear V for all historic data in dependence on the 
already used up lifetime, measured by the number A of flights. In the 
specific example, approximately 1300 flights have been made at the cur- 

20 rent point in time and the wear lies at approximately 55%. 



An extrapolation to the limit value G for the temperature increase T now 
takes place for the prediction of the future course of the wear V and thus 
for the estimate of the residual service life. There are naturally numerous 

25 possibilities for an extrapolation. The selection of a suitable extrapolation 
depends on the specific application. Consideration must in particular be 
taken of how conservative the estimate may or should be, that is how high 
specifications on the security of the forecast must be made. Accordingly, 
more or less pessimistic developments can be taken into account in the 

30 extrapolation and computed by the extrapolation. 
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In very sensitive applications, such as in the present case of an aeroplane 
engine, for example, one will work with more pessimistic forecasts in order 
to have sufficient safety reserves and to preclude a reaching of the wear 
5 limit prior to the end of the estimated residual service life. 

In the embodiment described here, a pessimistic development of the tem- 
perature increase can be derived over the number of flights, which are 
normally at least required to reach the limit value G = 40 K for the tem- 
10 perature increase T. It can be estimated, for example with reference to 

such wear curves Kl - K4, such as are shown in Fig. 1, what is the mini- 
mum number A of flights up to the reaching of the limit value. This mini- 
mum number is given in Fig. 1 by the point where the wear curve Kl 
reaches the limit value. 

15 

The minimum number of flights results at 2000 flights from the data 
underlying the code field KF in Fig. 3. 

The procedure is now as follows with respect to the extrapolation. Fig. 5 
20 shows the plane which is set up by the A axis and the T axis. As adready in 
Fig. 3, only one point (Ak, tk) is drawn. Let this be the most recent point, 
that is that point which corresponds to the present in which the estimate 
is made. From this point, as is shown by the chain dotted line E, an ex- 
trapolation is made in a linear manner to the limit value G = 40 K for the 
25 temperature increase T. The gradient with which the line E is "attached" to 
the point Tk, Ak), has been determined as follows. Starting from the fact 
that the minimum number of flights to reach the limit value G of the 
temperature increase T amounts, as mentioned above, to A = 2000 flights, 
the associated increase is determined which results when the points T = 0, 
30 A = 0 (new state) and T = 40 L, A = 2000 flights are connected to one 
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another by a straight line. The gradient resulting from this is then reduced 
again for safety reasons; in the present case it was halved, that is, it was 
assumed that the limit value is already reached after half the number of 
flights. The extrapolation resulting from this is the line E shown in Fig. 5. 

5 

In the next step, the line E is shown in the code field KF, whereby it gen- 
erally becomes a curved line. If the code field KF is now again projected - 
analogously to the representation in Fig. 4 - onto the plane set up by the 
V axis and by the A axis, the representation shown in Fig. 6 results. In 

10 addition to the historic values (which are naturally identical to those in 

Fig. 4), the projection E' of the extrapolation can now also be seen which is 
shown in a chain dotted manner. The projection E' of the extrapolation 
reaches the wear limit V = 1 after a number of flights which is designated 
by AG as the end value. The estimate of the residual service life now re- 

15 suits by forming the difference between the end value AG and the current 
value for the number of flights, that is the latest of the historic values. In 
Fig. 6, this is the last point of the continuous curve. 

This method for the estimating of the residual service life is now repeated 
20 constantly or at pre-settable intervals during the service life of the engine. 
The prediction horizon is substantially defined by the choice of the gradi- 
ent. The choice of a reasonable gradient has to take place in an applica- 
tion specific manner. 

25 For better understanding, Fig. 7 illustrates the estimating of the residual 
service life for seven stages in the life cycle of an engine. In each of the 
applications in Fig. 7 the plane is shown, in each case in the same man- 
ner as in Fig. 4 and in Fig. 6, which is set up by the V axis and the A axis 
and which results by projection of the code field. The historic data, that is 

30 those data which are based on data obtained by a technical measurement, 
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are each shown as continuous lines; the extrapolations which underlie the 
estimate are shown in a chain dotted manner. The transition from the 
continuous to the chain dotted line therefore gives the respective present 
in which the estimate took place. The estimated residual service life RL 
5 can be seen in number of flights in each case to the right next to the 

representations. In the topmost representation, in which no historic data 
are yet present, the residual service life is substantially defined by the 
chosen gradient. 

) 

10 In the second to fifth figures, seen from the top, a residual service life of 
approximately 500 flights is always forecast. This means that the aero- 
plane operator receives the information after, for example, 650 flights 
(third figure) that at least 500 flights are still possible up to the next ser- 
vice. After approximately 1000 flights (fourth representation), the operator 

15 receives the information that at least 480 flights are still possible. After 

somewhat more than 1300 flights (fifth representation), the operator again 
receives the information that 480 flights are still possible. These forecasts 
in the second to fifth representations are also substantially defined by the 
prediction horizon and thus by the choice of the gradient. 

) 20 

Only in the sixth representation (sixth estimate) does the end of the life 
cycle emerge so that the operator can now initiate the audit planning. Due 
to the selected gradient for the extrapolation, the emerging end of the life 
cycle can be recognised approximately 500 flights before reaching the wear 
25 limit. If this is related to the whole life cycle of approximately 2000 flights, 
substantial free room for action of approximately 25% of the life cycle 
remains for the operator. 

The prediction horizon can admittedly be expanded or extended by selec- 
30 tion of a smaller gradient for the extrapolation. This would, however, not 
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result in a much larger benefit in practice, but would increase the risk of 
too optimistic an estimate of the residual service life. 

It becomes possible by the estimating of the residual service life of the 
5 aeroplane engine which is described here and which takes the quantified 
wear into account to use the potential of the engines substantially more 
efficiently, without compromises in the operating safety being required. 
The total service or audit planning can be optimised, whereby an eco- 
j nomically much more favourable operation is made possible. The method 

10 in accordance with the invention can thus in particular be used advanta- 
geously for the service planning on aeroplane engines. This also makes a 
much more efficient planning of the service work possible at a plurality of 
aeroplanes. The method in accordance with the invention consequently 
makes an extremely high performance fleet management of a whole fleet of 
15 aeroplanes possible. 

The type of extrapolation described above is naturally only to be under- 
stood as an example. Different kinds of extrapolation, e.g. non-linear 
extrapolations or extrapolations based on qualitatively known develop- 
) 20 ments can also be used. The specific choice of an extrapolation matched to 
the respective application does not present the person skilled in the art 
with any problems. 

Even if, in the specifically described preferred embodiment for the deter- 
25 mination of the code field, a linguistic fuzzy model is used for the process- 
ing of a-priori know-how, the invention is not restricted to such modelling 
processes. 
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It is also not necessarily the case that the code field must be established 
by means of a-priori know-how. Other methods are also possible to deter- 
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mine a code field which gives the relationship between the characteristic 
parameter, the time parameter and the wear. For instance, depending on 
the application, for example ab-initio computations can be carried out, or 
interpretation calculations, dimensioning calculations, physical modelling 
processes, data supported modelling processes, system behaviour compu- 
tations. Furthermore, it is possible to describe the code field in the form of 
polynomials, look-up tables, multi-layer perceptrons (neuronal networks), 
radial basis functions, Singleton and Takadi-Sugeno fuzzy models as well 
as Hinging hyperplanes. 

It is naturally also possible to use more than one characteristic parameter 
sensitive to the wear, whereby the code field is defined in a higher dimen- 
sional space. 



